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FIGURES 

1  An  example  of  the  stress  versus  strain  response  for  the  composites  of  table  1  1 1 

for  most  temperatures  and  strain  rates.  The  dashed  line  indicates  schematically 

the  decrease  of  the  stress  at  the  lowest  temperatures  for  some  of  the  composites. 

2  Log  of  the  compressive  strength  versus  log  of  Young’s  Modulus  for  PAX  2A.  1 2 

The  straight  lines  were  drawn  in  with  the  indicated  slopes. 

3  Log  of  the  compressive  strength  versus  log  of  Young's  Modulus  for  LX-14.  1 3 

The  straight  lines  were  drawn  in  with  the  indicated  slopes. 

4  Strain  at  the  maximum  stress  versus  the  reciprocal  of  the  square  root  of  the  14 

modulus  for  PAX  2A.  The  straight  line  was  fitted  to  the  data  points  for  1/(Eo)0  6 

values  from  0  to  0.04  with  R2=0.878. 

5  Strain  at  the  maximum  stress  versus  the  reciprocal  of  the  square  root  of  the  modulus  1 5 

for  LX-14.  The  straight  line  was  fitted  to  the  data  points  lying  along  the  line  with 
R2=0.948. 

6  Stress  versus  strain  curves  normalized  to  the  maximum  stress  for  LX-14  at  65°C  1 6 

for  several  strain  rates. 

TABLES 

1  Composition  of  particulate  polymer  composites  17 
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INTRODUCTION 


This  work  was  initiated  to  survey  the  mechanical  properties  of  a  group  of  particulate  polymer 
composite  explosives  (refs.  1  through  4).  The  polymer  composites  are  made  up  of  polymer  binders 
(with  plastizer  in  most  cases)  and  80%  to  95%  organic  polycrystalline  nonpolymer  particulate  explo¬ 
sives  (table  1).  The  general  approach  of  the  work  presented  was  to  determine  the  stress-strain 
(stress  versus  strain)  properties  of  these  composites  as  a  function  of  temperature  and  strain  rate. 
The  condition  of  the  samples  after  deformation  was  also  noted;  i.e.,  whether  there  was  evidence  of 
plastic  deformation,  cracking,  and/or  fracture.  All  of  this  information  was  then  used  in  developing 
insight  into  the  nature  of  the  failure  processes. 


EXPERIMENTAL 

Stress-strain  data  in  compression  were  obtained  using  an  MTS  servo-hydraulic  system 
operated  at  constant  strain  rates  of  0.001  to  30/s  (ref.  5).  Samples  were  in  the  form  of  right  circular 
cylinders  34  in.  to  1  in.  in  length  and  34  in.  to  %  in.  in  diameter  and  the  end  faces  of  the  samples  were 
coated  with  a  lubricant,  e.g.  graphite,  to  minimize  frictional  effects  between  the  sample  and  the  load¬ 
ing  platens.  Samples  were  conditioned  at  temperatures  between  -60°  and  75°C  for  at  least  2  hrs 
before  measurements  and  were  then  compressed  along  the  cylinder  axis  to  obtain  engineering 
stress  and  strain.  One  to  five  samples  were  measured  at  each  temperature  and  strain  rate.  How¬ 
ever,  the  curves  of  figures  1  and  6  and  the  points  of  figures  2  through  5  each  represent  one  sample. 

Samples  of  the  polymer  composite  explosives  were  prepared  either  by  pressing  to  size  or  by 
pressing  into  large  billets  and  machining  to  size  (refs.  1 ,  2,  and  6).  Precautions  were  taken  to  insure 
that  the  cylinder  end  faces  were  flat  and  parallel.  The  densities  of  all  samples  were  measured  and 
results  are  presented  only  for  samples  having  densities  in  a  narrow  range  close  to  the  maximum 
theoretical  (zero  porosity)  density. 

In  table  1,  the  composition  of  the  explosive  composites  considered  are  given.  Glass  transition 
temperatures,  Tg,  are  given  where  known. 


RESULTS 

For  uniaxial  compression  of  the  materials  of  the  table,  the  stress  initially  increases  linearly  with 
increasing  strain  then  curves  over  and  passes  through  a  maximum  stress  with  further  increases  in 
strain  as  shown  in  figure  1 .  The  stress  either  decreases  continuously  for  additional  increases  in 
strain  beyond  the  maximum  in  the  strain-softening  region  (fig.  6),  or  in  some  cases,  at  the  lowest 
temperatures  the  stress  decreases  abruptly  to  near  zero  for  strains  at  or  just  beyond  the  maximum 
as  indicated  schematically  in  figure  1.  In  some  cases,  there  is  non-linearity  at  the  very  beginning  of 
the  stress-strain  curve  due  to  sample  end  conditions  (not  shown  in  figure  1).  Failure  is  considered 
here  to  initiate  at  the  point  given  by  of  and  s,  in  figure  1 ,  where  the  stress-strain  curve  deviates  from 
linearity.  Stresses  up  to  about  this  point  do  not  cause  permanent  changes,  while  stresses  beyond 
this  point  do  cause  permanent  changes  (ref.  7).  However,  this  point  is  often  difficult  to  locate 
precisely  because  of  noise  and  curvature.  Therefore,  the  maximum  stress,  the  compressive 
strength,  is  used  here  as  a  measure  of  the  failure  stress,  and  the  strain  at  the  maximum  stress  is 
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used  as  a  measure  of  the  failure  strain.  For  the  data  of  figure  1 ,  this  is  approximately  equivalent  to 
using  a  0.4%  strain  offset  as  a  definition  of  failure.  Three  quantities  taken  from  the  initial  part  of  the 
stress-strain  curves  are  of  interest: 

•  The  initial  slope,  which  is  taken  as  a  measure  of  the  undamaged  modulus,  E0 

•  The  maximum  compressive  stress  (the  compressive  strength),  am,  taken  as  a 
measure  of  the  failure  stress 

•  The  strain  at  the  maximum  stress,  sm,  taken  as  a  measure  of  the  failure  strain 

(fig- 1). 

After  strains  much  greater  than  em,  the  samples  appear  to  be  plastically  deformed  and/or  macro- 
scopically  cracked,  and  in  some  cases  fractured.  Fractured,  as  used  here,  means  broken  into 
separate  pieces.  In  those  cases  at  lower  temperatures  where  the  stress  has  dropped  abruptly  to 
near  zero  at  or  near  the  maximum  (fig.  1),  the  samples  are  always  fractured,  and  in  some  cases 
(e.g.,  PAX  2  and  PAX  2A)  they  are  fragmented  into  many  pieces. 

In  figures  2  and  3,  the  log  of  the  compressive  strength,  om,  is  given  versus  the  log  of  the 
modulus,  Eo,  for  PAX  2A  and  for  LX-14,  respectively.  Data  is  given  for  temperatures  between  -45°C 
to  65°C  and  for  strain  rates  between  0.001  Is  and  1.0/s.  For  the  smaller  values  of  the  modulus,  the 
compressive  strength  increases  linearly  with  the  modulus  while  for  the  larger  values,  the  strength 
increases  as  the  square  root  of  the  modulus.  Linear  plots  (not  shown)  of  am  versus  E0  (refs.  2  and  3) 
and  am  versus  E01/2  for  the  appropriate  regions  of  figures  2  and  3,  give  straight  lines  that  pass 
through  the  origin.  These  results  suggest  two  separate  failure  modes  or  branches:  a  failure  mode 
with  the  failure  strength,  am,  proportional  to  E0  for  smaller  values  of  E0,  designated  here  as  the  linear 
failure  mode  and  a  failure  mode  with  cm  proportional  to  E01/2  for  larger  values  of  E0,  designated  here 
as  the  square-root  failure  mode.  This  conclusion  is  supported  by  the  observations  that  the 
compressive  strength  is  exponen-tially  dependent  on  1/temperature  (1/K)  with  two  activation 
energies:  one  for  the  region  of  the  linear  failure  mode  and  another  for  the  region  of  the  square-root 
failure  mode  (ref.  8).  The  latter  activation  energy  is  one  half  of  the  former.  The  modulus  is  also 
exponentially  dependent  on  1 /temperature  (1/K),  but  with  the  same  activation  energy  throughout 
both  regions  and  with  a  value  equal  in  magni-tude  to  the  value  for  the  strength  in  the  linear  mode 
(ref.  8).  The  temperature  and  strain  rate  depen-dencies  of  the  modulus  are  thought  to  be  primarily 
associated  with  the  viscoelastic  properties  of  the  binder. 

For  a  given  value  of  E0,  it  appears  that  the  failure  mode  (linear  or  square  root)  is  determined 
by  the  mode  that  requires  the  lower  value  of  failure  stress,  i.e.,  the  compressive  strength.  Therefore, 
for  values  of  log  E0  less  than  about  2.65  for  PAX  2A  (fig.  2),  the  linear  mode  requires  the  lower 
stress  while  for  values  of  log  E0  greater  than  2.65  the  square  root  mode  requires  the  lower  stress. 
Similar  considerations  apply  to  the  data  for  LX-14  in  figure  3.  The  linear  mode  was  observed  for  all 
of  the  composites  of  the  table,  but  the  square  root  mode  was  observed  only  for  PAX  2,  PAX  2A,  LX- 
14,  and  to  a  limited  extent  for  9404  and  Comp  A3  Type  II.  All  of  the  available  data  for  9501  lies  in 
the  linear  range,  while  the  available  data  for  9502  for  the  larger  modulus  range  does  not  conform  to 
either  a  linear  or  a  square  route  relationship.  N9  is  anomalous  in  that  the  stress-strain  curve  before 
the  maximum  (fig.  1),  indicates  two  approximately  linear  regions  and  so  two  slopes.  The  slope  at 
the  larger  strain  is  greater  than  the  initial  slope.  The  peak  stress  is  proportional  to  the  initial  slope 
and  so  indicates  the  linear  failure  mode.  Several  other  composites  and  TNT  also  exhibit  the  linear 
failure  mode  (ref.  3).  The  square  root  branch  may  not  have  been  observed  to  date  in  some  of  the 


composites  because  the  stress  and  the  moduli  required  may  be  greater  than  those  encountered 
experimentally.  Lower  temperatures  and/or  higher  strain  rates  may  then  be  necessary  to  observe 
this  branch  in  these  composites. 

From  the  stress-strain  curve  of  figure  1 ,  a  geometric  relationship  between  the  three  quantities 
cm,  E0,  and  em  can  be  obtained  as 

<^m  =  Eo  /("I  "*■  a)  (1) 

where  cm  (1  +  a)  is  the  stress  at  which  a  straight  line  through  the  initial  portion  of  the  stress-strain 
curve,  the  slope  of  which  defines  E0,  intersects  a  constant  strain  line  at  sm  (ref.  2).  The  parameter  a 
is  a  measure  of  the  shape  of  the  stress-strain  curve  between  the  point  where  it  deviates  from  a 
straight  line  and  the  point  of  maximum  stress.  In  all  cases  where  a  linear  relationship  between  cm 
and  E0  is  observed,  sm  is  found  to  be  approximately  constant  (refs.  1  through  4)  as  illustrated  in 
figure  4  for  PAX  2A  and  in  figure  5  for  LX-14  for  the  smaller  values  of  E0  (larger  values  of  1/E01/2). 
While  the  parameter  a  does  change  somewhat  with  temperature  and  strain  rate,  the  magnitude  of 
the  changes  in  a  are  such  that  equation  1  is  satisfied  with  an  approximately  constant  sm  for  all  the 
materials  of  the  table  and  in  the  temperature  and  strain  rate  range  such  that  cm  is  proportional  to  E0. 
For  PAX  2A,  the  parameter  a  has  a  maximum  as  a  function  of  temperature  in  the  low  temperature 
range  and  the  temperature  of  the  maximum  increases  with  strain  rate.  This  may  contribute  to  the 
apparent  scatter  of  the  data  of  figure  4  in  the  low  temperature  (large  E0)  range.  The  data  for  LX-14 
and  the  other  materials  are  not  sufficiently  extensive  to  determine  if  similar  maxima  in  the  parameter 
a  occur. 

Although  em  is  approximately  constant  (compared  to  the  changes  in  cm  and  E0)  with  changes  in 
temperature  and  strain  rate  in  the  ranges  for  the  linear  failure  mode,  there  is  a  tendency  for  this 
strain  for  some  of  the  materials  to  increase  slightly  with  decreasing  temperature  and  increasing 
strain  rate. 

Substituting  am  =  K  E01/2  into  equation  1  for  the  square  route  branch  and  rearranging  gives 

£m  =  K  (1  +  a)/E01/2  (2) 

This  equation  is  satisfied  for  the  appropriate  values  of  E  and  for  those  composites  which  exhibit  the 
square  route  branch.  This  is  also  illustrated  for  PAX  2A  in  figure  4  and  for  LX-14  in  figure  5  for  the 
larger  values  of  E  (smaller  values  of  1/E1/2). 


DISCUSSION 

Studies  of  the  uniaxial  stress-strain  behavior  following  prior  uniaxial  deformation  (and  so 
damage)  indicate  that  for  maximum  prior  strains  in  the  initial  linear  region  of  the  stress-strain  curve 
(fig.  1),  a  change  in  the  modulus  is  not  detected.  However,  for  maximum  prior  strains  exceeding  this 
linear  range,  the  modulus  decreases  continuously  with  increasing  prior  strain  and  so  amount  of 
damage  (ref.  7  and  9).  These  results  indicate  that  the  stress-strain  response  of  figure  1  and  well  into 
the  softening  region  (fig.  6)  can  be  attributed  to  a  damage  modulus  that  decreases  continuously  with 
increasing  strain  beyond  the  linear  range  (refs.  7  and  9).  Therefore,  the  uniaxial  stress-strain 
response  can  be  expressed  as 
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a  =  E  s  =  (E0  -  AE)  e  =  e/Y  =  s/(Y0  +  AY) 


(3) 


where  E  is  Young's  modulus,  Y  is  the  compliance,  E0  and  Y0  are  the  values  before  damage  due  to 
loading,  arid  AE  and  AY  are  the  changes  due  to  damage  introduced  by  loading.  To  first  order  AE/E0 
=  AY/Y0.  AY  is  assumed  to  be  a  single-valued  monotonicaliy  increasing  function  of  e.  Y0  and  AY 
may  be  functions  of  temperature,  strain  rate,  initial  crack  and  dislocation  parameters,  density,  etc. 
However,  all  damages  (e.g.,  changes  in  crack  and  dislocation  parameters)  is  assumed  to  be  ex¬ 
pressed  through  the  dependence  of  e  on  these  changes.  Dienes  (ref.  10)  and  Dienes  and  Riley 
(ref.  11)  used  equation  3  with  several  models  for  the  dependence  of  AY  on  crack  parameters  in 
fitting  the  author's  data  for  9501 . 

By  differentiating  a  of  equation  3  with  respect  to  e,  and  setting  the  resulting  differential  equal  to 
zero,  the  condition  for  a  maximum  in  the  curve  of  o  versus  s  is  found  to  be 

sm  =  (Y0  +  AYm)/AYm'  (4) 

where  the  subscript  m  denotes  the  values  at  the  maximum  and  the '  symbol  indicates  differentiation 
with  respect  to  s,  e.g.,  AYm'  =  [d(AY)/de]m.  By  substitution  of  equation  4  into  3,  the  stress  at  the 
maximum  is 


om  =  1/AYm'  (5) 

From  equation  4,  AYm'  must  be  non-zero  and  positive.  Thus,  AYm  must  be  an  increasing  function  of 
s  as  assumed. 

To  satisfy  the  linear  branch  of  the  am  versus  E0  curve,  it  is  then  necessary  for  AYm'  to  be 
expressed  as 


AYm'  =  Y0f(sm)'  =  f(em)7E0  (6) 

with  f(em)'  independent  of  or  insensitive  to  temperature  and  strain  rate.  It  is  assumed  here  that  the 
initial  linear  slope  of  the  stress-strain  curve  is  given  by  E0.  Substitution  of  equation  6  into  4  yields 


sm  =  (Y0  +  AYm)/Y0f(em)' 


(7) 


If 


AYm  =  Y0f(em)  (8) 

em  is  independent  of  or  insensitive  to  temperature  and  strain  rate  in  agreement  with  the  experiment. 
Thus,  the  relationships  for  am  and  em  are 

am=Eo/f(sm)'  (9a) 


and 
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em  =  (1+f(Sm))/f(em)' 


(9b) 


where  f(sm)  and  f(em)'  are  independent  of  or  insensitive  to  temperature  and  strain  rate. 

If  f(s)  is  valid  at  all  strains  and  is  expressed  as  a  power  series  in  s,  the  condition  that  there  is  a 
maximum  in  the  a  versus  s  curve  is  that  the  highest  power  of  s  is  greater  than  unity.  The  total 
compliance  from  equations  3  and  8  is  then  given  by 

Y  =  Y0(1+f(s))  (10) 


so  that 


(1+f(e»  (11) 

is  a  damage  function  for  the  compliance. 

In  summary,  by  expressing  the  total  compliance  as  the  sum  of  the  undamaged  compliance 
and  the  change  in  compliance  due  to  damage,  and  by  comparing  the  resulting  equations  for  am  and 
sm  with  the  experimental  results,  it  is  concluded  that  the  change  in  compliance  is  proportional  to  the 
undamaged  compliance  (eq.  8)  and  further  that  the  damage  function  is  primarily  a  function  of  the 
total  strain  and  is  independent  of  or  insensitive  to  temperature  and  strain  rate.  Other  conclusions 
may  be  possible. 

The  linear  branch  of  the  am  versus  E0  curve  and  the  attendant  constant  £m  have  been  dis¬ 
cussed  elsewhere  in  terms  of  either  a  constant  strain  as  the  criterion  for  failure  or  in  terms  of  an 
unspecified  mechanism  requiring  om  to  be  proportional  to  E0  (refs.  2  and  3).  From  equation  1 ,  a 
constant  em  results  in  am  being  proportional  to  E0  and  conversely  if  am  is  proportional  to  E0,  sm  must 
be  constant,  all  if  the  parameter  a  is  constant.  This  damage  approach  provides  a  rationale  for  a 
constant  strain  criteria  for  failure  with  changes  in  temperature  and  strain  rate  (eq.  9b)  based  on  a 
damage  function  that  is  independent  of  or  insensitive  to  these  two  parameters  in  the  initial  stages  of 
damage  (eq.  11).  In  addition,  the  linear  relationship  between  om  and  E0  is  obtained  (eq.  9a).  There¬ 
fore,  this  damage  approach  satisfies  the  experimental  criteria.  This  damage  mechanism  also,  of 
course,  provides  the  rationale  for  why  am  has  the  same  temperature  and  strain  rate  dependencies 
and  so  activation  energy  as  E0  in  the  linear  failure  mode  range  (ref.  8).  As  noted,  the  temperature 
and  strain  rate  dependencies  of  E0  are  thought  to  be  due  primarily  to  the  viscoelastic  properties  of 
the  binder. 

Results  at  25°C  and  a  low  strain  rate  also  indicate  that  if  the  strain  does  not  exceed  a  strain  of 
about  sf  (fig.  1),  a  change  in  the  modulus  within  the  precision  of  the  measurements  is  not  detected. 
However,  if  sf  is  exceeded,  the  modulus  is  decreased  as  noted  (ref.  7  and  9).  Therefore,  the  dam¬ 
age  function  (eq.  11)  may,  at  least  for  practical  purposes,  have  a  damage  threshold  ef.  The  thres¬ 
hold  stress,  Of,  then  varies  linearly  with  the  modulus  as  temperature  and  strain  rate  vary.  Further 
work  in  this  area  is  desirable. 

While  additional  work  is  necessary  to  establish  the  nature  of  the  damage,  the  conclusions  that 
the  change  in  compliance  is  proportional  to  the  undamaged  compliance  and  that  the  damage  func¬ 
tion  is  independent  of  or  insensitive  to  temperature  and  strain  rate  places  significant  restrictions  on 
the  possibilities.  Thermally-activated  slow-crack  growth  is  not  a  possibility  because  the  function  f(s) 
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(eq.  8)  is  not  a  function  of  temperature  and  strain  rate  (refs.  5  and  8).  Brittle  fracture  of  particulate 
particles,  perhaps  involving  particle  interactions,  may  under  some  circumstances  satisfy  the  neces¬ 
sary  conditions.  Evidence  for  decreased  particulate  particle  sizes  in  several  of  the  composites  of  the 
table  after  damage  by  uniaxial  compression  were  found  by  small  angle  x-ray  and  neutron  scattering 
(ref.  12).  While  the  author  is  unaware  of  mechanical  property  measurements  for  any  of  the  com¬ 
posites  of  the  table  as  a  function  of  particulate  particle  size,  measurements  were  made  of  pressed 
HMX,  the  particulate  for  most  of  the  composites  of  table  1,  for  two  different  average  particle  sizes 
(ref.  13).  Samples  containing  only  coarse  HMX  gave  moduli  and  compressive  strengths  higher  by 
about  30%  compared  to  samples  containing  a  bimodal  mixture  of  coarse  and  fine  particles  sizes  in 
the  ratio  3/1 .  Both  sample  groups  were  pressed  to  approximately  the  same  density.  In  contrast, 
measurements  of  samples  of  an  inert  simulant  of  9501  containing  sugar  in  place  of  HMX,  and  with 
the  same  two  particle  size  distributions  as  for  the  HMX  discussed  previously,  did  not  show  differ¬ 
ences  in  the  moduli  or  the  compressive  strengths  (ref.  14).  Measurements  of  the  moduli  and  the 
particle  sizes  for  the  composites  of  the  table  as  a  function  of  deformation  or  other  variables  are 
desirable  to  further  resolve  this  matter. 

A  normalized  stress-strain  relationship  can  be  obtained  from  equations  3,  5,  and  8  as 

a/am  =  Ef(8m)'/(1+f(s))  (12) 

Therefore,  to  the  extent  that  the  equations  used  to  obtain  equation  12  are  valid  at  strains  other  than 
sm,  this  normalized  stress-strain  relationship  is  independent  of  or  insensitive  to  temperature  and 
strain  rate  if  f(e)  and  f(em)'  are  independent  of  or  insensitive  to  temperature  and  strain  rate.  An 
examination  of  the  data  indicates  that  the  normalized  stress-strain  curves  are  relatively  insensitive  to 
temperature  and  strain  rate  in  the  initial  linear  region  and  in  the  vicinity  of  the  maximum,  but  are 
more  sensitive  to  these  variables  in  the  softening  region.  This  is  illustrated  in  figure  6  by  normalized 
stress-strain  curves  for  LX-14  at  65°C  for  several  strain  rates.  The  differences  in  the  softening 
region  indicate  a  rather  strong  strain  rate  effect  while  the  slopes  (moduli)  in  the  initial  region  and  the 
strains  at  the  maxima  indicate  a  minimal  strain  rate  dependence.  The  strain  rate  dependence  in  the 
softening  region  decreases  with  decreasing  temperature  and  is  minimal  at  the  lowest  temperatures 
(not  shown).  Some  of  the  differences  in  the  initial  region  before  the  maxima  are  due  to  deviations  of 
the  sample  end  surfaces  from  flatness  and  parallelism.  In  general,  similar  results  were  obtained  for 
other  composites  of  the  table.  However,  it  is  important  to  note  that  the  samples  are  not  damaged 
uniformly  along  the  sample  length  in  the  softening  region.  After  loading  into  the  softening  region, 
the  samples  are  often  barreled,  i.e.,  the  permanent  radial  expansion  is  less  at  the  sample  ends  than 
in  the  mid-section.  This  effect  is  due  to  the  frictional  constraints  on  the  radial  expansion  of  the 
sample  ends  by  the  compression  platens.  This  occurred  even  though  lubrication  was  used  on  the 
sample  ends  to  minimize  this  friction.  Therefore,  the  observed  stress-strain  curves  are  averages 
over  ranges  of  damage.  The  softening  region  of  the  stress-strain  curves  will  be  treated  in  more 
detail  elsewhere  (ref.  15). 

Dienes  and  Riley  (ref.  11)  have  used  equation  3  with  Y  expressed  in  the  form  of  equation  10 
and  f(s)  =  k  e3,  as  obtained  from  a  very  simple  crack  propagation  model,  to  fit  this  author’s  data  for 
9501  for  a  range  of  temperatures  and  strain  rates.  Y0  (E0)  and  k  were  varied  to  give  the  best  fits  to 
the  stress-strain  curves  for  strains  between  zero  and  0.04,  the  latter  being  well  into  the  softening 
region  of  the  curves.  The  fits  to  the  data  with  k  and  E0  dependent  on  temperature  and  strain  rate  are 
surprisingly  good  considering  the  simplicity  of  the  model.  However,  the  model  predicts  a  much 
stronger  dependence  of  k  on  strain  rate  than  obtained  by  fitting  to  the  data  and  the  fits  to  the  initial 
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portions  of  the  stress-strain  curves  of  immediate  concern  here  are  in  several  cases  not  as  good  as 
might  be  desired.  The  model  predicts  sm  to  be  dependent  on  k  and  so  temperature  and  strain  rate. 

In  addition,  the  model  predicts  am  to  be  proportional  to  E0  as  observed,  but  also  predicts  am  to  be 
dependent  on  k.  In  particular,  this  dependence  of  <rm  on  k  and  so  temperature  and  strain  rate,  in 
addition  to  the  dependence  of  cm  on  temperature  and  strain  rate  through  E0,  is  contrary  to  observa¬ 
tions  (figs.  2  and  3).  The  temperature  dependence  of  k  is  not  discussed  by  Dienes  and  Riley  (ref. 
11).  As  noted,  the  normalized  curves  are  insensitive  to  temperature  and  strain  rate  for  the  initial 
parts  of  the  curves  but  are  temperature  and  strain  rate  dependent  in  the  softening  region  (fig.  6).  A 
more  complex  model  may  be  necessary  to  fit  the  data  in  the  initial  regions  and  in  the  softening 
region.  Dienes  and  Riley  also  fit  some  of  this  author's  9501  data  to  a  more  complex  crack  model 
with  realistic  crack  parameters,  but  the  fits  are  not  as  good  as  the  fits  with  the  simpler  model  (ref. 

11).  The  goal  of  this  latter  fitting  was  to  deduce  information  on  the  microstructure,  e.g.,  crack  size 
and  number  density.  These  are  used  in  a  hot-spot  analysis  of  ignition  associated  with  mechanical 
deformation.  Unfortunately,  these  authors  used  a  form  of  the  data  for  all  of  their  curve  fitting  that 
was  not  corrected  for  measuring  system  compliance.  This  correction  influences  the  data  primarily  at 
lower  temperatures  (larger  stresses). 

Aidun  has  developed  a  phenomenological  viscoelastic  model  for  9501  based  on  maxwell 
elements  (ref.  15).  A  damage  function  was  deduced  to  account  for  the  difference  between  this 
author's  9501  data  at  25°C  and  the  viscoelastic  model.  Reasonably  good  fits  were  obtained  over 
essentially  the  entire  stress-strain  curves  by  using  a  damage  function  which  is  strain  rate  indepen¬ 
dent  in  the  initial  parts  of  the  curves  up  to  about  the  maxima  (fig.  1),  but  which  is  strain  rate  depen¬ 
dent  in  the  softening  region  at  greater  strains.  This  strain  rate  dependence  of  the  damage  function 
is  consistent  with  this  discussion  of  the  parameter  k  used  by  Dienes  and  Riley  (ref.  1 1).  Tempera¬ 
ture  variations  were  not  considered  by  Aidun. 

In  addition,  Gozonas  has  developed  a  nonlinear  viscoelastic  model  with  damage  due  to  slow 
crack  growth  which  he  has  fitted  to  M30  propellant  data  as  a  function  of  temperature  and  strain  rate 
for  strains  in  the  vicinity  of  and  greater  than  sm  (refs.  16  and  17).  M30  has  similar  mechanical 
properties  to  the  materials  considered  here  (ref.  3).  The  agreement  between  the  model  and  the  data 
is  reasonably  good  for  this  range  of  strains  and  does  include  the  observed  constant  failure  strain,  em 
(refs  3  and  16).  However,  the  model  was  not  fitted  to  the  initial  linear  part  of  the  stress-strain  curve 
defining  E0  and  does  not  appear  to  predict  this  part  of  the  curve  very  well.  Therefore,  the  model 
does  not  predict  the  observed  linear  relationship  between  am  and  E0  (ref.  3). 

The  approach  used  here  for  the  linear  mode,  based  on  a  damaged  compliance,  is  not  readily 
adaptable  to  the  square  root  failure  branch  of  the  am  versus  E0  curve.  There  are,  however,  other 
possible  reasons  for  the  square  root  relationship  between  cm  and  E0  not  directly  involving  damage 
functions,  but  instead  failure  criteria.  Consider  first  the  total  work  per  unit  volume,  Wv,  performed  in 
uniaxial  compression.  Wv  can  be  expressed  as 

Wv  =  Jcds  (13) 

Equation  13  can  be  written  as 

Wv  =  1/EoJ  ada  =  (af)2/2E0  (14) 

when  the  integration  is  over  the  initial  linear  portion  of  the  stress-strain  curve  from  zero  stress  to  a 
failure  stress  af  taken  at  the  point  where  the  curve  (fig.  1)  deviates  from  linearity.  Then  from 
equation  14 
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cf=(E0)1/2(2Wv)1/2 


(15) 


is  obtained.  For  an  elastic  material  Wv  is  equal  to  the  total  elastic  energy  per  unit  volume,  Uv. 
However,  because  the  materials  in  question  were  found  to  have  viscoelastic  properties  (ref.  9),  the 
elastic  energy  will  be  a  fraction  r  of  Wv.  Therefore,  equation  15  can  be  written  as 

<jf  =  (E0)1/2  (2Uv/r)1/2  (16) 

With  <jm  taken  as  a  measure  of  af,  equation  16  gives  the  observed  relationship  between  om  and  E0  if 
(Uv/r)  is  constant  as  temperature  and  strain  rate  are  varied  over  the  limits  of  the  square  root  region  of 
the  Gm  versus  E0  curve  (figs.  2  and  3).  For  elastic  materials  (r=1),  equation  16  gives  the  relationship 
between  the  failure  stress  and  E0  for  a  total  maximum  elastic  strain  energy  criterion  for  failure  (ref. 
18).  A  similar  relationship  is  obtained  for  the  elastic  energy  of  distortion  criterion  and  the  equivalent 
maximum  octahedral  shearing  stress  (von  Mises)  criterion  for  failure  (ref.  1 8).  But,  because  r  is  not 
constant  for  viscoelastic  materials,  it  appears  that  these  failure  criteria  are  not  applicable  to  the 
materials  in  question.  However,  before  concluding  that  these  failure  theories  are  not  in  fact  applic¬ 
able  because  of  variations  in  r,  it  is  necessary  to  determine  experimentally  the  magnitude  of  the 
changes  in  r  over  the  range  of  temperature  and  strain  rate  of  interest. 

Additional  measurements  are  necessary  to  distinguish  between  the  several  failure  criteria 
mentioned  here.  It  should  be  possible  to  distinguish  between  a  constant  total  strain  energy  criteria 
and  a  constant  strain  energy  of  distortion  criteria  by  comparing  the  results  of  torsion  tests  with  the 
uniaxial  results  (ref.  18).  The  total  strain  energy  criteria  is,  however,  unlikely  because  failure  has  not 
been  observed  in  other  materials  for  hydrostatic  compression  to  total  strain  energies  significantly 
higher  than  the  value  for  uniaxial  failure  (ref.  18). 

Consider  next  the  Griffith  condition,  which  in  contrast,  is  related  to  the  decrease  of  the  total 
strain  energy  due  to  crack  growth  and  occurs  when  this  decrease  is  equal  to  the  increase  in  the 
surface  energy  also  due  to  crack  growth.  Frictional  effects  between  fracture  surfaces  are  neglected. 
The  Griffith  condition  may  be  stated  as 

=  [y  E0  K/c  ]1/2  (17) 

where  cg  is  the  stress  required  for  crack  growth,  c  is  the  crack  length,  y  is  the  surface  energy  per  unit 
area,  and  K  is  a  numerical  constant  (ref.  19).  While  equation  17  is  valid  for  a  particular  set  of 
conditions  which  differ  from  the  experimental  conditions  used  here,  the  relationship  between  the 
stress  and  the  modulus  should  be  generally  valid  since  it  is  obtained  from  the  relationship  for  the 
strain  energy.  Because  of  this  relationship,  the  slope  of  the  square  root  branch  (figs.  2  and  3)  is 
then  a  function  of  the  (maximum)  crack  length  c  for  this  mechanism  of  failure.  This  dependence  on 
the  maximum  crack  length  could  explain  the  large  scatter  of  the  experimental  points  of  figure  2  for 
this  branch  because  of  variations  of  maximum  crack  length  from  sample  to  sample.  Measurements 
of  failure  strength  as  a  function  of  crack  length  are  desirable  to  determine  if  the  square  root  branch  is 
indeed  due  to  the  Griffith  condition  for  crack  growth.  A  characteristic  of  this  crack  growth  (without 
crack  arrest)  is  abrupt  fracture  as  indicated  by  a  sharp  decrease  of  the  stress  to  zero  or  near  zero 
with  increasing  strain.  Although  abrupt  fracture  (fig.  1)  is  observed  at  the  lowest  temperatures  (large 
E0)  for  most  of  the  composites  exhibiting  the  square  root  branch,  it  is  not  observed  for  all  of  the 
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composites  having  this  branch  and  is  not  observed  for  many  of  the  points  lying  on  this  branch  (figs. 

2  and  3).  Thus,  if  the  square  root  branch  is  due  to  the  Griffith  condition  for  crack  growth,  it  must  be 
concluded  that  crack  arrest  is  also  important  in  many  cases  because  of  this  lack  of  abrupt  fracture. 
Damage  must  be  accumulated  by  this  crack  growth  and  arrest  and  in  some  cases  apparently  leads 
to  eventual  abrupt  fracture  at  the  lowest  temperatures. 

It  is,  therefore,  possible  to  interpret  the  initial  part  of  the  stress-strain  response  in  terms  of 
deformation-induced  damage  that  is  detected  as  increased  compliance.  The  incremental  compli¬ 
ance  (due  to  damage)  is  proportional  to  the  undamaged  compliance  and  the  damage  function  is 
determined  either  exclusively  or  primarily  by  the  total  strain  and  is  independent  of  or  insensitive  to 
temperature  and  strain  rate.  This  interpretation  applies  to  the  composites  of  table  1  and  some  other 
materials  (ref.  3),  all  of  which  give  linear  relationships  between  the  maximum  stress,  am,  and  the 
modulus,  E0  and  a  constant  or  almost  constant  sm.  In  addition,  as  just  considered,  some  of  the  com¬ 
posites  of  table  1  appear  to  exhibit  another  mechanism  which  produces  failure  at  larger  E0  (lower 
temperatures  and/or  higher  strain  rates)  at  significantly  lower  stresses  than  the  stresses  required  for 
this  damage  process.  In  these  cases,  om  varies  as  the  square  root  of  E0.  This  type  of  relationship 
can  be  associated  with  the  Griffith  conditions  for  crack  growth,  with  a  critical  total  strain  energy  for 
failure  or  with  a  critical  strain  energy  of  distortion  for  failure.  It,  therefore,  appears  that  for  some  of 
the  composites,  there  is  a  shift  from  a  strain-dependent  damage  mechanism  to  a  stress  dependent 
failure  process  as  the  temperature  is  lowered  and/or  the  strain  rate  is  increased.  The  results  also 
indicate  that  for  any  given  set  of  conditions  of  temperature  and  strain  rate  and  so  modulus  E0,  the 
failure  process  which  is  dominant  (i.e.,  with  am  proportional  to  either  E0  or  E01/2)  is  the  process  that 
requires  the  lower  stress.  Additional  work  is  required  to  further  establish  the  nature  of  these  failure 
processes. 


SUMMARY 

The  results  suggest  two  failure  modes  or  branches  in  the  temperature  and  strain  rate  ranges 
studied,  a  lower  modulus  mode  in  which  the  failure  strength  is  proportional  to  the  modulus  and  the 
failure  strain  is  (approximately)  constant,  and  a  higher-modulus  mode  in  which  the  failure  strength  is 
proportional  to  the  square  root  of  the  modulus  and  the  failure  strain  decreases  as  the  modulus 
increases.  The  lower  modulus  mode,  which  was  observed  for  all  of  the  composites  studied,  occurs 
at  higher  temperatures  and  lower  strain  rates  while  the  higher  modulus  mode,  which  was  observed 
only  for  some  of  these  composites,  occurs  at  lower  temperatures  and  higher  strain  rates.  The 
results  suggest  that  the  lower  modulus  mode  can  be  attributed  to  a  damage  compliance  with  the 
change  in  compliance  proportional  to  the  undamaged  compliance  and  a  damage  function  that  is 
strain  dependent  and  is  independent  of  or  insensitive  to  temperature  and  strain  rate.  The  results 
further  suggest  that  the  criterion  for  failure  for  the  higher  modulus  mode  is  either  a  constant  strain 
energy  or  a  critical  stress  for  crack  growth.  Additional  work  is  indicated. 
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Figure  1 

An  example  of  the  stress  versus  strain  response  for  the  composites  of  table  1  for  most  temperatures  and  strain  rates. 
The  dashed  line  indicates  schematically  the  decrease  of  the  stress  at  the  lowest  temperatures  for  some  of  the  composites. 
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Figure  2 

e  compressive  strength  versus  log  of  Young's  Modulus  for  PAX  2A. 
The  straight  lines  were  drawn  in  with  the  indicated  slopes. 
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Figure  3 

Log  of  the  compressive  strength  versus  log  of  Young's  Modulus  for  LX-14. 
The  straight  lines  were  drawn  in  with  the  indicated  slopes. 
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Figure  4 

Strain  at  the  maximum  stress  versus  the  reciprocal  of  the  square  root  of  the  modulus  for  PAX  2A 
The  straight  line  was  fitted  to  the  data  points  for  1/(E0)°5  values  from  0  to  0.04  with  R2=0.878. 
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Figure  5 

Strain  at  the  maximum  stress  versus  the  reciprocal  of  the  square  root  of  the  modulus  for  LX-14. 
The  straight  line  was  fitted  to  the  data  points  lying  along  the  line  with  R2=0.948. 
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Figure  6 

Stress  versus  strain  curves  normalized  to  the  maximum  stress  for  LX-14  at  65°C  for  several  strain  rates 


Table  1 

Composition  of  particulate  polymer  composites 


Binder 

Tg 

Name 

Particulate 

Polvmer 

Plastizer 

(°C) 

1  Pax  2 

HMX 

CAB 

BDNPA/F 

-37a 

80% 

8% 

12% 

II  Pax  2A 

HMX 

CAB 

BDNPA/F 

-37a 

85% 

6% 

9% 

III  9404 

HMX 

NC* 

CEF 

-Q 

CO 

94% 

2% 

3.84% 

IV  9501 

HMX 

ESTANE 

BDNPA/F 

-41  (B)° 

95% 

2.5% 

2.5% 

V  9502 

TATB 

KEL  F  800 

_Q 

o 

CO 

95% 

5% 

VI  LX-14 

HMX 

ESTANE  5702-F1 

-31  (B)b 

95.5% 

4.5% 

VII  Comp  A3 

RDX 

POLYETHYLENE 

91% 

9% 

VIII  N9 

HMX 

HYCAR  4454 

DOA 

-43d 

92% 

2% 

6% 

Nomenclature:  HMX  -  Cyclotetramethylene  tetranitramine. 

TATB  -  1 ,3, 5-triamino-2, 4, 6-trinitrobenzene. 

RDX  -  Cyclotrimethylene  trinitramine. 

NC  -  Nitrocellulose.  NG  -  Nitroglycerine. 

CAB  -  Cellulose  Acetate  Butyrate. 

BDNPA/F  -Bis(2,2-Dintropropyl)Acetal/Formal. 

CEF  -  Tris(Beta  Chloroethyl)  Phosphate. 

Estane  -  Polyurethane. 

KEL  F  800  -  Chlorotrifluoroethylene/vinylidine  flouride  copolymer. 
DOA  -  Dioctyl  adipate. 

B  -  Property  of  the  binder. 

*  Also  contains  0.1%  Diphenylamine 

a  Personal  Communication,  J.  Harris,  Picatinny  Arsenal,  N.J. 
b  Reference  20. 
c  Reference  21. 
d  Reference  22. 
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